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Abstract 

We analyze subsamples of Abell and ACO cluster catalogs, in order to study the 
spatial properties of the large scale matter distribution. The subsamples analyzed 
are estimated to be nearly complete and are the standard ones used in the corre- 
lation analysis by various authors. The statistical analysis of cluster correlations is 
performed without any assumptions on the nature of the distribution itself. The 
cluster samples show fractal correlations up to sample limits (~ Idh^Mpc) with 
fractal dimension D ~ 2, without any tendency towards homogenization. Our anal- 
ysis shows that the standard correlation methods are incorrect, since they give a 
finite correlation length for a distribution that does not possess one. In particular, 
the so called correlation length ro is shown to be simply a fraction of the sample 
size. Moreover we conclude that galaxies and clusters are two different representa- 
tions of the same self similar structure and that the correlations of clusters are the 
continuation of those of galaxies to larger scales. 



1 Introduction 



The determination of the matter distribution in space represents a crucial test 
for the initial hypothesis of the main cosmological theories [|]]. Both stan- 
dard Friedmann and Steady State models assume, beyond a certain scale Ao, 
the homogeneity of the matter distribution. This assumption can be tested, 
properly investigating the spatial distribution of galaxies and clusters. 

With respect to the galaxy catalogs, cluster surveys offer the possibility to 
study the large scale structure of the matter distribution in volumes much 
larger, reaching depths beyond z ~ 0.2 and extending all over the sky. Using 
clusters, one can trace matter distribution with a lower number of objects 
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with respect to the galaxies in the same volume. For example, in the northern 
hemisphere we know pa 10 6 galaxies up to z pa 0.15 that correspond to pa 500 
rich clusters. However the problem of cluster catalogs is their incompleteness, 
since clusters are identified as density enhancements in galaxy angular surveys 
and their distance is usually determined through the redshift measurement of 
one or two galaxy member. It is clear in fact that only the measurements of 
every clusters (and member galaxies) would allow to construct truly volume 
limited samples. Moreover, as we show in the following, there is a strong 
arbitrariness in the definition of a cluster, such an arbitrariness can be avoid 
using the methods of modern statistical mechanics. 

Cluster distribution shows large inhomogeneities and huge voids. Tully (0, 
||), investigating the spatial distribution of Abell and ACO clusters up to 
300/i -1 Mpc, observes the presence of structures on a scale of 0.1c, lying in 
the plane of the Local Supercluster. Galaxies clump in clusters and clusters 
in superclusters with extensions comparable to the largest scales of current 
samples. Catalogs of superclusters, compiled by various authors (Q, 0, ||, 
0, H) show the presence of correlated structures with extension ranging from 
few Mpc to 150h~ 1 Mpc. Complementary to the presence of large structures 
is the existence of voids. Several authors (f5|, [|], Q, flOfl , ||]) have investi- 
gated the shape and the dimension of voids determined by rich clusters and 
superclusters. The dimension of voids is usually defined as diameter of empty 
sphere containing no clusters. As for clusters, the definition of voids implies 
a certain degree of arbitrariness, since one can consider only spherical voids 
or ellipsoidal voids and so on. In particular Einasto et al. (||) find that in 
their cluster sample (up to z ~ 0.1), the median radius of voids is 50/i _1 Mpc. 
However, the observed voids have elongated shapes and in one dimension they 
can exceed this value. The authors observe in fact a giant void of 300/i _1 Mpc 
of length. A more detailed study on a single void is performed by [|ll|], which 
have investigated the distribution of galaxies in and around the closest one, 
the Northern Local Void. The authors found that the dimension of the voids 
depends on the objects which have been used to defined them: voids defined 
by clusters are larger than galaxy defined voids, which are larger than faint 
galaxy defined voids. There is, then, a hierarchy of voids. However voids seem 
to be real and not due to observational incompleteness. Search of extremely 
faint dwarf galaxies in voids has given, up to now, negative results ( |T2|, see 
refs in [[[]]]). Moreover one observes an increase of the void size with sample 
size, that may correspond to a a self similarity in the distribution of voids. 
Even in this case we show that the modern methods of statistical mechan- 
ics allow a study of void distribution does not depend on any assumption of 
arbitrariness. 

The observation of large scale structures like superclusters and voids, extend- 
ing up to the limits of the samples, clearly make questionable the existence 
of an average density: the average density may be not a well defined intrinsic 
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property of the system and the assumption of homogeneity for the distribution 
must be checked and not assumed. Usually, the presence of inhomogeneities 
on the scale of sample analyzed is interpreted as incompleteness of the samples 
themselves. However, up to now, better and more extensive observations have 
confirmed the reality of these structures, making the observed agglomeration 
sharper and not filling the voids. 

An argument usually proposed to reconcile the observed inhomogeneities to 
the supposed large scale homogeneity, is that the structures are indeed real, 
but their amplitudes are small; in other words the density contrast 5N/N ,i.e. 
the ratio of the density fluctuations with respect the average density, tends to 
zero at the limits of the sample ( J0|, [p6| , jl5| ). One expects, therefore, that 



going a bit further, the homogeneity would be eventually reached. A part from 
the fact that this expectation has been systematically disproved, the argument 
is conceptually wrong. As reported in 5N/N tends to zero at the limits of 
the sample for any system and can not be interpreted as tendency toward 
homogenization. This is because the average is computed within the sample 
volume. 

The standard measure of the irregularities in the space distribution is the 
dimensionless autocorrelation function (CF) £(r). This is the density-density 
correlation normalized by the average density of the system. If the average 
density is not an intrinsic property of the system, unique and independent, 
the results are spurious ([ffH], [0])- The basic point is that the use of the 
standard CF assumes implicitly that the homogeneity is well reached within 
the sample size, and consequently it can not test whether it happens or not. 

For both galaxy and cluster distributions, £(r) function is found to be a power 
law: 



£(r) « (r /r) 7 (1) 



both the exponent and the amplitude r are uncertain. The exponent 7 gener- 
ally appears consistent with the value 1.8 observed for galaxy autocorrelation 
function ( |19|], [20|| , [ pl|1 , p2f , p3| ). The quantity r is the so called "correlation 



length" and is defined as £(ro) = 1. 

In general high values for r (20 — 25h~ 1 Mpc) are obtained for samples that 
contains more rich Abell/ACO clusters (f2|, f§, g|, |6|, |7|, @, [||]), 



although there is a debate whether these high values might be produced by 
systematic biases present in the Abell/ACO catalog. The proposed corrections 
for these effects produce a lower r of « Uh^Mpc ( J30J, 0, @, 0). 
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In conclusion, for the Abell/ACO samples analyzed so far: 

Uh^Mpc <r < 25h~ 1 Mpc (2) 



According to the standard analysis, r should be the correlation length for the 
distribution and its physical meaning is the following: below this value, clus- 
ters should be strongly correlated, while correlations should became small and 
eventually negligible at scale of the order of twice this length. At this distance 
the distribution should became rather smooth and homogeneous. The obser- 
vational evidences go instead exactly in the opposite direction, since larger 
samples show larger correlated structures. The result of the standard statis- 
tical analysis, i.e. a little correlation length, seems in contradiction to the 
observational evidences. Moreover clusters are found to be more clustered 
than galaxies, for which the standard analysis obtains a correlation length tq 
hhr x Mpc ( J35[| , |55fl). This mismatch between galaxy and cluster correla- 



tions is a puzzling feature of the usual analysis; clusters, in fact, are made 
by galaxies and many of these are included in the galaxy catalogs for which 
ro ~ 5h~ l Mpc is derived. It is therefore necessary to assume that cluster 
galaxies have fundamental differences with respect these that are not part of 
a cluster. This concept has given rise to the so called richness clustering re- 



lation (0, p7| , [pSfl , p9|); according to this, objects with different mass or 
morphology would segregate from each other and produce different correlation 
properties; in other words, the correlation length of a given class of objects i 
is rQi m 0.4 < n > i , where < n >j is relative mean spatial density. This 
hypothesis has been applied to explain the increase of r$ going from APM and 
X — ray less rich clusters to the Abell more rich clusters. Szalay and Schramm 
40| have interpreted the scaling of r as signature of scale invariance prop- 



erties of the distribution. However, as pointed out by [0 (hereafter CP92), 
this is not correct since for a self-similar distribution tq is proportional to R s , 
where R s is the effective sample radius (see sec. 5) and not to < n > -1 / 3 . In 
effect, statistical analysis of sparser objects requires larger sampling volumes. 
As the following analysis will show, cluster distribution is fractal up to the 
sample size; then we can reinterpret the scaling of tq as an increase of R s . 



The criticisms of the standard analysis have been raised by Pietronero and 
collaborators ( |I6|| , |T7| , CP92). In CP92, the authors analyze, with methods of 
modern statistical mechanics, the CfAl catalog of galaxies ( f4T|) and a sample 
of clusters from Abell catalog ( [0). This analysis, differently from the £(r), 
has no a priori assumptions on the nature of the distribution and is commonly 
used in modern statistical mechanics, where one has to face non-analytical 
distributions. The results are that the galaxy and cluster distributions are 
not homogeneous within the sample limits, but on the contrary, fractal up to 
the sample limits for statistical analysis: these are ~ 20h~ 1 Mpc for galaxy 
and to ~ 60 — 80/i _1 Mpc for cluster in the samples analyzed. The analysis 
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is performed to greater distances for clusters simply because the volume of 
cluster sample is larger than the galaxy one. In other words, there is not 
any scale of homogeneity and correlation properties of clusters are just the 
continuation of the galaxy correlations to larger scales. The mismatch between 
the cluster and galaxy correlation is due to a mathematical inconsistency (the 
use of and it is automatically eliminated using the correct statistical 

analysis. 

Recently such an analysis has been performed also on several other galaxy 
surveys, namely Perseus Pisces redshift catalog, LEDA, ESP: these results 
confirm and extend those of CP92. Galaxy distribution in Perseus Pisces cat- 
alog is not homogeneous, but fractal with the fractal dimension D pa 2 up to 
30k- 1 Mpc ( @, H, @, U); in LEDA database, fractal up to lhOh^Mpc 



with dimension D pa 2 ( ||47|| ). ESP catalog requires a slight different anal- 
ysis, which implies a lower degree of statistical validity. With these remarks, 
ESP catalog shows a fractal behaviour with the same dimension (D pa 2) up 
800 — 900/i _1 Mpc ( 0). Putting all these results together, one obtains power 
law (i.e. fractal) correlations for galaxy distribution, extending from few Mpc 
up to 800 — §Wh~ 1 Mpc PS| . Consequence of fractal correlations is that r , 



derived from £(r), is not the correlation length, but simply a fraction of the 
effective sample radius R s . According to this, scales in all samples with R s . 

In these paper we present results of the statistical analysis on various sam- 
ples of Abell/ACO clusters. In section 2 we give a description of the samples 
analysed. In section 3 we discuss the samples completeness and the behaviour 
of the density versus radial distance. In section 4 we briefly introduce the 
basic properties of fractal distributions. In section 5 we present the results 
of the correlation function T(r) analysis. In section 6 we study the standard 
£(r) function and we analyse the sample depth dependence of the correlation 
length ro and finally we summarize the results and draw our conclusions. 



2 Description of the data and subsamples 



The study of the spatial distribution of clusters requires a complete sample of 
clusters over a large volume. There are many available extensive catalogs of 
clusters from which one may try to extract complete subsamples. 

Here we have analyzed subsample extracted from Abell catalog and ACO 
catalog. From Abell catalog we have analyzed the Postman sample ( 0); this is 
a sample of 351 clusters with the tenth ranked galaxy magnitude (mio)< 16.5. 
The sample includes all such clusters which lie north of 5 = — 27°30'. The 
typical redshift of a cluster with mxo < 16.5 is z ~ 0.09. To this redshift, 
incompleteness in the Abell catalog is considered insignificant ( 0). About 
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half of clusters have the redshifts based on at least three independent galaxy 
spectra, 25% have redshift determined from two galaxy spectra and the last 
25% from a single spectrum ( 0). 

Postman et al. (1992) [[F] have defined a so called statistical subsample, which 
consist of 208 clusters with z < 0.08 and \b\ > 30°. This subsample has been 
shown to be minimally affected by selection biases. 

Bahcall & Soneira 1983 sample ( |J, BS83 hereafter): this sample includes 104 
clusters of distance class D < 4 (z < 1), richness class R > 1 and located at 
high galactic latitude > 30°) 



From ACO catalog |27| we have selected all the cluster with measured redshift 
in according to the following constraints: mio < 16.4,6 < —20°, 5 < — 17° ( 
50| , H49| ). The sample contains 139 clusters up to the limiting distance of 



930h~ 1 Mpc. 



3 Sample completeness 



The completeness of the samples is usually estimated from the density of 
the clusters as function of the redshift. In Fig.l, it is shown the density in 
the samples analysed. The density n(r) shows large fluctuations followed by 
a power law decay as ~ r~ 3 . Computation of the density in bins produce 
larger fluctuations, because this is a differential quantity more noisy than the 
integral density n(r). Usually the fluctuating behaviour of the density up to 
z < 0.1 is interpreted as a flat one, i.e. the density of the sample is considered 
to be constant and the distribution homogeneous; consequently the sample 
is considered complete and homogeneous up to this distance and incomplete 
beyond. The decay r -3 is, in fact, due incompleteness of the samples, since the 
number of clusters is nearly constant, but the volume sampled increases. We 
note that up to the beginning of the incompleteness region, all the samples 
contain few clusters; the sparsest is the BS83 sample: the northern galactic 
part contains 53 clusters up to 230h~ 1 Mpc and the southern part 24. The 
Postman sample |7| contains 120 clusters up to z ~ 0.08 in the northern part 
and 88 in the southern one. The ACO sample 91 up to the 150/i -1 Mpc. The 
basic point is that n(r) is not an averaged quantity, i.e. is the density computed 
from the vertex. In this situation at small scales, there are no clusters and the 
density is zero; going a little bit further, one starts to see some clusters, and 
the density shows large fluctuations because the statistics is low. At larger 
distance, the fluctuations reduce for the increasing of the statistics, but the 
average behaviour is apparent only sampling a quite large range of scale. In 
a homogeneous distribution, the amplitude of fluctuations reduce with scale, 
while in a fractal these have the same amplitude in a log scale, because the 
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Fig. 1. The radial density n(r) from the vertex for a) Abell catalog Postman sample, 
north (diamonds) and south (filled circles) b) Abell catalog BS83 sample north 
(diamonds) and south (filled circles) c) ACO catalog sample (empty circles). 

structure is self-similar and is made by fluctuations (see sec. 4). Hence, to 
recover the mean behaviour, one has to sample a quite large range of scale. 

In Fig.l we see that the range of scale, over which we can measure n(r), is 
quite small (~ 50h~ 1 Mpc) and in our opinion is not enough to recover the 
average behaviour of the density. 

The study the properties of cluster spatial distribution requires an almost 
complete sample; then it is necessary to exclude or correct the incompleteness 
region. One way to obtain this in the standard approach is to assume an 
homogeneous distribution of clusters up to the sample limits. The observed 
distribution in the redshift space is then weighted with a selection function 
p(z), that is the ratio between the observed counts of clusters in the volume 
dV(z) at redshift z and those expected from an homogeneous distribution. In 
our analysis, we want to avoid any assumption on the distribution itself and 
for this reason we will limit our analysis up to a depth corresponding to the 
beginning of the incompleteness region, without correcting by means of any 
selection function. 

Another selection effect must be taken in account, regarding to the observed 
depletion of the surface mean density of clusters at low galactic latitude (to 



< 30°). This is probably due to obscuration and confusion with high-density 



regions of stars of our galaxy ( ||26|| ). As in the case of redshift incompleteness, 
one way to overcome this incompleteness is to weigh the observed distribution 
with a latitude selection function P(b), that is the ratio between observed sur- 
face cluster density at latitude b and the expected from an homogeneous one. 
The normalization of this selection function is arbitrary, because it depends 
from the real density. Another way is to limit the sample to high galactic lat- 
itude region. We will adopt this standard procedure (i.e. \b\ > 30° for Abell 
catalog and b < —20° for ACO), that has no assumptions with only the in- 
convenient to slightly limit the sample. 

For the Postman sample the distance of completeness is estimated to be 
roughly 230h~ 1 Mpc corresponding to z < 0.08, for the BS83 sample the same 
distance, while for ACO sample is ~ 150h~ 1 Mpc. All these limits have been 
estimated from the behaviour of the density versus the distance (Fig. 1). 



4 Properties of fractal distributions 



In this section we mention the essential properties of fractal structures because 
they will be necessary for the correct interpretation of the statistical analysis. 
However in no way these properties are assumed or used in the analysis it- 
self. A fractal consists of a system in which more and more structures appear 
at smaller and smaller scales and the structures at small scales are similar 
to the ones at large scales. The self similarity of these structures is then in- 
compatible with analyticity. Standard mathematical tools based on analytical 
functions can not characterize these distributions. The first quantitative de- 
scription of these forms is the metric dimension. One way to determine it, is 
the mass-length method. Starting from an point occupied by an object, we 
count how many objects N ("mass") are present within a volume of linear 
size r ("length") ( |J): 

N(r) =B-r D (3) 



D is the fractal dimension and characterizes in a quantitative way how the 
system fills the space. The prefactor B depends to the lower cut-offs of the 
distribution; these are related to the smallest scale above which the system is 
self-similar and below which the self similarity is no more satisfied. In general 
we can write: 
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where r* is this smallest scale and N* is the number of object up to r*. For 
a deterministic fractal this relation is exact, while for a stochastic one it is 
satisfied in an average sense. 

Eq.(^|) corresponds to a smooth convolution of real N(r), that is a very fluc- 
tuating function; a fractal is, in fact, characterized by large fluctuations and 
clustering at all scales. We stress that eq. (§) is valid in general; for an homo- 
geneous distribution, for example, one has D = 3. 

From eq. (|3]) , we can compute the average density < n > for a sample of radius 
R s which contains a portion of the structure with dimension D. The sample 
volume is assumed to be a sphere ( V(R S ) = (4/3)-7ri?f) and therefore 

< n >- BSA - LbbtW (5) 

< >_ V(R S ) ~ 4n a [b) 



If the distribution is homogeneous, D = 3 and the average density is constant 
and independent from the sample volume; for a fractal, the average density 
depends explicitly on the sample size R s and it is not a meaningful quantity. 
It is a decreasing function of the sample size and < n >^ for R s — > oo. It 
is important to note that eq. @ holds from every point of the system, when 
considered as the origin. This feature is related to the non-analyticity of the 
distribution. In a fractal distribution every observer is equivalent to any other 
one, i.e. it holds the property of local isotropy around any observer ( Sylos 
Labini 1994 |52[). We can define the conditional density from an occupied 
point i as: 

r Wi = S(r)-^ = Rbt-»-» (6) 



where S(r) is the area of a spherical shell of radius r. T(r)i is then the 
density at distance r from the i — th point in a shell of thickness dr. As for 
eq.(||), eq.® corresponds to a smooth convolution of fluctuating quantity, 
i.e. the conditional density from one element of the distribution. Usually the 
exponent that defines the decay of the conditional density (3-D) is called the 
codimension and it corresponds to the exponent 7 of the galaxy distribution. 



5 The conditional average density 



The correlation function suitable to study homogeneous and inhomogeneous 
distribution is described by ( |TB|; CP92) 

G(r) =< n{r + fi)n{fi) >j~ r~ 7 (7) 
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where the exponent 7 = 3 — D (in 3- dimensional space) (Eq. |6|) and the index 
% of the average means that this is performed on all the occupied points of 
the system. If the sample is homogeneous,/} = 3, G(r) ~< n > 2 and then 
is constant; if the sample has correlations on all scales, it is fractal, D < 3, 
7 > and G(r) is a power law. For a more complete discussion we refer the 
reader to CP92. We can normalize the G(r) to the size of the sample under 
analysis and define, following CP92: 

r /x = < nif+r^njfj) > { = G(r) 
<n> <n> 



where < n > is the average density of the sample. This normalization does 
not introduce any bias even if the average density is sample-depth dependent, 
as in the case of fractal distributions (Eq.||), because it represents only an 
overall normalizing factor. The T(r) (Eq. ||) can be computed by the following 
expression 



1 N 1 



r(r) = — - — 7r-r— / n(fi + r f )dr' 
i=i r . 



1 N D 

= ±y T (r) l = —Br 3 - D (9) 
N ~^ 47r y ! 

where iV is the number of objects in the sample. This is called the conditional 
average density (CP92). T(r) is the average of T(r)j and hence it is a smooth 
function away from the lower and upper cutoffs of the distribution (r* and the 
dimension of the sample). From eq.(^j) we see that T(r) is independent from 
the sample size, depending only by the intrinsic quantities of the distribution 
(B and D). It is also very useful to use the average density 

r 

r*(r) = / A7rr' 2 T(r')dr' (10) 

47rr d J 




This function produce an artificial smoothing of T(r) function, but it correctly 
reproduces global properties (CP92). 

As we said, T(r) (and r*(r)) is the average density computed in spherical 
shell. In such a way we eliminate from the statistics the points for which a 
sphere of radius r is not fully included within the sample boundaries. This 
prescription allows us to avoid any weighting scheme, i.e. any assumption in 
the treatment of the boundaries conditions. Of course in doing this, we have 
a smaller number of points and we stop our analysis at a smaller depth than 
that of other authors (0,0, j2"9"fl) In fact, we have to limit our analysis 



10 





N 




< n > 




J. Lg 


' 






sr 


h 3 Mpc~ 3 


(h^Mpc) 


(h^Mpc) 


(h^Mpc) 


BS83 sub north 


53 


3.1 


4.2 • 10~ 6 


230 


70 


27 ±2 


Postman s.s. north 


120 


3.1 


9.5 • 10~ 5 


230 


70 


26 ±2 


Postman s.s. south 


58 


1 i 
i. / 


J..O • 1U 




Do 




ACO sub 


91 


2.2 


3.7- 10~ 5 


150 


50 


17 ± 1 



Table 1 

features of the various samples analyzed; column I: the catalog from which the 
sample has been extracted; column II: number of clusters in the sample; column 
III: solid angle of the sample; column IV: density of the sample; column V: sample 
extension; column VI: radius of the maximum sphere fully included; column VII: 
the correlation length 

to an effective depth R s that is of the order of the radius of the maximum 
sphere fully contained in the sample volume (CP92; see also [48|). We have 
studied the behavior of T(r) and T*(r) in the samples of Table 1. The results 
are shown in Fig. 2. In Fig2a we have reported F*(r) for the north and south 
galactic part of Postman statistical subsample. 

The south galactic part has a smaller solid angle with respect the north one 
and by consequence a smaller R s . A well denned power law behavior is detected 
up to the sample limit without any tendency towards homogenization. The 
codimension is, with good accuracy 
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3-D w 1.0 ±0.1 



so that D p& 2.0 ± 0.1 up ~ 70/i~ 1 Mpc for the north part and up to w 
50h~ l Mpc for the south. The result is in good agreement with those from 
BS83 sample (Fig. 2b). The north galactic part is fractal with dimension D w 
1.7 ±0.2 up to r w 70h~ 1 Mpc. We have not reported the analysis for the BS83 
south galactic part, because it gives a very noisy result, due to the very poor 
statistics of the sample. The lower statistics of the BS83 sample with respect 
Postman one is also the reason for the little difference in the estimation of 
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r(Mpc) 

Fig. 2. The behaviour of the conditional density T(r) and average conditional den- 
sity T*(r) for a) Abell catalog , Postman sample respectively north (diamonds) and 
south galactic parts (filled circles) b) Abell catalog ,BS83 sample north galactic part 
(diamonds), c) ACO catalog sample (filled circles). T(r), T*(r) are power law up to 
the radius of the greatest sphere fully included in the sample. The reference line has 
slope 7 1 



the fractal dimension D of the two samples. We note that the codimension 7 
for these samples has a lower value with respect the standard determination; 
we refer the reader to next section for a explanation of this effect. In Fig. 2c 
we show the results from ACO sample. The T(r) is a power law with expo- 
nent D ~ 2.0 ± 0.1 up to bOh^Mpc. The fluctuations at r < lOh^Mpc 
are due to the fact that these distances are of the order of the minimum 
average distance between clusters (ss Yhhr l Mpc for Postman and BS83 and 
~ 8h~ 1 Mpc for ACO sample), so that we can then interpret them as a effect 
of poor statistic. In the correct regime, at distances > 10/i -1 Mpc, where the 
samples are statistically significative, the slope of T(r) is almost the same for 
the Postman sample and for the ACO and, with a slight difference because 
the poor statistics, for the BS83 sample. All the samples investigated have, 
then, consistent statistical properties, i.e. they show a clear behaviour of the 
correlation function. A power law behaviour corresponds to long range corre- 
lations, which cannot be produced by random incompleteness in the sample, 
if they are not present. Long range correlations (fractal) can be only destroyed 
by incompleteness, but not produced by it. 
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Hence the samples show well defined statistical properties, i.e. they are statisti- 
cally fair samples; the results of the analysis is that they are not homogeneous 
samples, but, on the contrary, fractal. 



6 £(r) analysis 



CP92 clarify some crucial points of the standard CF analysis, and in particular 
they discuss the meaning of the so-called "correlation length" r found with 



the standard approach ( |21|) and defined by the relation: 



*(ro) = 1 (12) 
where 

(r) = <n(r>(rl + r)>, _ i 
<n> 2 

is the two point correlation function used in the standard analysis. As we said, 
if the average density is not a well defined intrinsic property of the system, 
the analysis with £(r) gives spurious results. In particular, if the conditional 
density T(r) (or T*(r)) is a power law, the system is fractal and the average 
density is simply related to the sample size. In this case, if one derives a 
correlation length r , comparing the average density-density correlation to 
the square of average density, one obtains simply a fraction of the effective 
sample radius R s . In other words, if a system has a correlation function T(r) 
that is a power law, the system is self-similar and it has no reference values 
(like the average density) with respect to which one can define what is big or 
small. Hence, computing a correlation length with respect the average density 
is meaningless. As pointed out in Baryshev et al. 1994 |J, this is true also 
for a system that has a characteristic length, as for example a system with 
correlation function T(r) = T e~ r ^ r °. In this case, the characteristic length 
r is related to the behaviour of the function and not the pref actor T ; the 
characteristic length is not defined by the condition T(ro) = 1. 

Following CP92, the expression of the £(r) in the case of fractal distribution, 
is: 

e(r) = ((3- 7 )/3)(r/ J R s )^-l (14) 



where R s (the effective sample radius) is the depth of the spherical volume 
where one computes the average density from Eq. (||). From Eq. (14) it follows 
that 
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i.) the so-called correlation length r$ (defined as £( r o) = 1) is a linear function 
of the sample size R s 

ro = ((3- 7 )/6)^i2 a (15) 



and hence it is a quantity without any correlation meaning but it is simply 
related to the sample size. 

ii. ) the amplitude of the £(r) is: 

A{R S ) = {{3- 1 )/3)W S (16) 

iii. ) £(r) is a power law only for 

((3-7)/3)(r/i2.)"^»l (17) 



hence for r < r^. for larger distances there is a clear deviation from a power 
law behavior due to the definition of £(r). However usually £(r) is fitted with 
a power law in the range 0.5ro ^ r < 2ro. The result is one obtains a greater 
exponent than 3-D. This is the reason why the usual estimation of the 
exponent 7 of £(r) is ~ 1.7, different from 7 » 1, corresponding to -D ~ 2, 
that we found by mean of the T(r) analysis. Moreover we point out that there 
is another spurious effect in the calculation of the £(r). The £(r) is estimated 
by means of the expression 

where Ndd(r) is the number of data pairs at separation r and N pp (r) is the 
number of pairs at separation r for a random distribution with the same den- 
sity and geometry of the survey. Eq.(18) corresponds to compute the density- 
density correlation not only in spheres fully contained in the sample volume, 
but also in portions of sphere. This produce an artificial homogenization, for 
distances r greater than the radius of the greatest sphere fully contained in 
the sample. For these separations, in fact, the density -density correlation is, 
of course, computed only in portions of sphere and this corresponds to assume 
that the density, found in the portion, is the same in the all solid angle. To 
avoid this effect, we have computed the £(r) only up to R s . We have studied 
the £(r) in our samples of clusters. In fig. 3 we have reported the £(r) for the 
Postman sample (north and south) and for the ACO sample. We have fitted 
the experimental points with the functional form of eq. (|i~4"|) . We found that 
7 « 1 for the all samples analyzed. The corresponding r$ are reported in the 
Table 1; for the Postman sample north we found r ~ 26 ± 2 and for BS83 



14 



10 1 



10~' 



10 J 



10' 




10 

r(Mpc) 



Fig. 3. The £(r) function for the a) Postman sample (north) (filled circles) b) Post- 
man sample (south) (filled circles) and ACO (diamonds). The reference line is the 
functional form of eq.(14) with 7 ~ 1 

sample tq ~ 27 ± 2 These two samples have the same as we expect, noting 
that they have the same R s . The values found are in agreement with eq.(15). 
The same conclusions hold for the other two samples: Postman south and 
ACO. Both of them have R s 50h~ 1 Mpc) and the same r (~ 17/i _1 Mpc), 
according to eq.(15). In conclusion r is simply a fraction of the sample size 
R s , without any physical meaning for the properties of the samples. 



7 Conclusion 



We have analyzed the properties of the various Abell/ACO cluster samples. 
The statistical analysis, performed without any a priori assumptions, shows 
that Abell samples are fractal with fractal dimension D rs 2 for Postman 
sample and for BS83 sample up to pa 70h~ 1 Mpc and for ACO sample up to 
pa 50/i _1 Mpc. The different limiting distance of the analysis in the various 
samples corresponds to the radius of the greatest sphere fully included in the 
sample, that is different in the various samples. This limitation allows us to 
avoid any weighting schemes or assumptions in the analysis. The result of the 
statistical analysis is that no tendency towards homogenization is detected 
within the sample limits. The so-called correlation length r derived from the 
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£(r) analysis, is simply directly proportional to the sample size R s and then, 
it is meaningless with respect the correlation properties of the system. 

This result is in agreement with the correlation analysis, performed measuring 
the conditional average density, T(r), that is a power law (fractal) within 
the sample limits. Moreover our results on cluster samples are in remarkable 
agreement with the same analysis performed on various galaxy redshift surveys 
(CfAl, Perseus Pisces, LEDA, ESP) as one expects for a fractal distribution 
of galaxy and clusters. The mismatch between galaxy and cluster correlation 
is just due to the mathematical inconsistency of the use of £(r) and the correct 
analysis, in terms of T(r) and r*(r), shows that correlations of clusters are 
just the continuation at larger scales of galaxy correlations. We conclude that 
galaxies and clusters are two different representations of the same self similar 
structure. 

Galaxy clusters extend the correlations of galaxies to deeper depth. In this re- 
spect cluster distribution represents a coarse grained representation of galax- 
ies; it is the same self-similar distribution, but sampled with a large scale 
resolution, i.e. considering a cluster of galaxies as a single object, without 
distinguish the structure in it. Hence, we can study the clusters distribution 
simply performing a coarse graining on galaxy distribution. Usually clusters 
are identified with some criteria, that are different according to different ob- 
servers. On the contrary, in this way, we can make the analysis independently 
from the definition of cluster, supercluster etc., but one has to know the com- 
plete galaxy survey over which performing the coarse graining procedure. The 
same considerations hold, of course, for the void distribution: the void distri- 
bution is, in fact, just the complement of the matter one. In conclusion, our 
methods have the advantage to be independent from the nature of the dis- 
tribution considered, and they gives a quantitative way to detect self similar 
properties, whether they exist, at different scales. 
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